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1 .0 SUMMARY 

■ 

V This document provides a set of preliminary entry guidance and autopilot soft- 

ware formulations for use in the Mission Control Center (MCC) entry processor. 
These software formulations meet all level B requirements as specified in refer 
ence 1 . 


2.0 INTRODUCTION 

This internal note presents the level C software formulations requirements for 
the entry guidance and the simplified autopilot that will be used by the MCC 
entry processor. Revision 2 incorporates the modifications required to function- 
ally simulate optional TAEM targeting capability (OTT). Implementation of this 
logic in the MCC must be coordinated with flight software OTT implementation (CR 
12848H) and MCC TAEM guidance OTT (Revision 2 to MCC Level C Formulation Require- 
ments for TAEM Guidance and Flight Control; JSC IN 76-FM-87). The entry guidance 
logic is based on the Orbiter avionics entry guidance software as described in 
reference 2. This MCC requirements document contains a definition of coordinate 
systems (3.2), a list of parameter definitions for the software formulations 

(3.3) , a description of the entry guidance detailed formulation requirements 

(3.4) , a description of the detailed autopilot formulation requirements (3*5), 
a description of the targeting routine (3-6)» and a set of formulation flow 
charts (appendixes A through C) . 


3.0 SOFTWARE FORMULATION REQUIREMENTS 


3.1 INTRODUCTION 

The entry guidance system is the source of the bank angle and angle-of-attack 
commands used to control the entry trajectory. The entry guidance can be called 
by one of two means, either in the normal mode to generate a complete entry 
trajectory or as a part of the iterative targeting mode in the entry target gen- 
eration (ETG) subphase of the deorbit processor. This logic will be discussed 
in section 3*4. 

The autopilot generates the Orbiter attitude response to the entry guidance 
commands. This is accomplished by means of a simple phase plane in the bank 
angle and angle-of-attack axis. The sideslip angle (6) is always assumed to be 
zero. This logic is discussed in section 3*5. 

The targeting logic generates the range and heading information to the targeted 
runway and is used by the entry guidance and the deorbit processor. This logic 
is described in section 3.6* 
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3.2 COORDINATE SYSTEMS 

Two basic coordinate systems are assumed by the software formulations described 
in this document. The state vector is assumed by the targeting routine (EGRT) 
to be in the Greenwich true-of-date system as defined in figure 3.2-1. The 
runway coordinate system is defined in figure 3.2-2 and is used by EGRT. The 
bank command and the angle-of-attack command generated by the entpr ^idance and 
executed by the autopilot are attitudes defined with respect to the Earth rela- 
tive velocity vector. 

3.3 PARAMETER LISTS AND DEFINITIONS 

A complete list of all parameters used in the entry guidance, autopilot, and 
targeting routine is presented with appropriate definitions. Tables 3.3-1 and 

3.3- 2 present the input data for the entry guidance. Table 3.3-3 presents a 
list of the internal parameters for entry guidance. Table 3.3-4, 3.3-5, and 

3.3- 6 present the same data for the autopilot and tables 3*3-7j 3*3 » 

3.3- 9 present the data for the targeting routine. 


3.14 ENTRY GUIDANCE FORMULATION 


3.4.1 Requirements Overview 

The entry guidance is the source for bank and angle-of-attack commands that are 
used to control the entry trajectory. The entry guidance can be called by 
either of two modes. If the entry guidance flag (EGFLG) is set to zero 
a normal entry guidance function is exercised, which will duplicate the Orbiter 
avionics entry guidance function. This mode will be used to simulate an entry 
trajectory and will also be used in the final iteration mode of the entry target 
generator (ETGl processor. The second mode (EGFLG=2) is used in the entry 
processor to simulate an entry based on a "canned" drag profile for the ETG 
targeting processor. 

In the normal mode (EGFLG=0), the entry guidance controls the entry trajectory 
by bank angle modulation while using a preselected angle-of-attack profile that 
is a function of relative velocity. Range predictions are based on solutions to 
the equation of motion for a specified entry drag-velocity profile. The drag- 
velocity profile, the shape of which is specified by the mission constants 
table, consists of quadratic, pseudoequilibrium glide, linear, and constant drag 
segments. Downrange errors are nulled by changing the magnitude of the Dank 
angle, and crossrange errors are nulled by bank angle reversals that limit the 
crossrange error within a converging deadband. This mode begins at 400 000 feet 
and ends at TAEM interface, relative velocity = VTAEM. 


In the bTG canned mode (EGFLG=2), the entry guidance controls the entry 
trajectory to a predefined drag-velocity profile. The entry simulations begin 
at 400 000 feet and are terminated after the pullout maneuver has been completed 
and the flight profile stabilized on the drag velocity profile at exactly 23 000 
ft/sec relative velocity (VETG). 
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In the normal mode, the entry guidance consists of five major phases: the 
preentry phase, the temperature control phase, the equilibrium glide pha^, 
constant drag phase, and the transition phase, as shown in figure 


the 


The preentry phase maintains the vehicle attitude until a p rede ter^ned tote 1 
load factor level is reached (ASTART). During this phase, the vehicle is 
maintained in a three-axis attitude hold mode. This attitude is specified by 
the flight controller as MM 304 <1,, a, and 0. This phase is terminated at 5.66 

ft/sec2 total acceleration level and the temperature control phase is entered. 


The temperature control phase is designed to control the entry trajectonr 
through pullout to a temperature profile consistent with the overall P*'°: 

file shape optimization and entry ranging requirements. This phase consists o. 
two quadratic drag-velocity segments that are used to optimize the entry pro- 
file? Range predictions are based on analytic solution of the equations of mo- 
tion for these two segments. This phase is terminated at the velocity, VBl, 
specified in the mission constants table, and the guidance is transferred to the 
equilibrium glide phase at this point. 


The equilibrium glide phase produces an equilibrium glide type trajectory, con- 
sistent with the ranging solution, until the equilibrium glide trajecto^ nnint" 
sects the constant drag trajectory required to reach the target. At this point, 
control is transferred to the constant drag phase, and range predictions are 
based on a constant drag profile until the transition phase te entered. At a 
specified velocity, VTRAN, (mission constants table), control is transferred to 
the transition phase. The transition phase is based on a linear drag profile, 
as a function of energy, which is required to null the range 

tion from the entry guidance to the TAEM guidance occurs at a velocity, VTAbM, 
specified in the mission constants table. Presently, this transfer point is at 
an Earth-relative velocity of 2500 ft/sec. 


The entry guidance generates bank angle and angle-of-attack commands to be used 
bv the autopilot. The bank angle commands are designed to converge the actual 
drag acceleration level to the reference drag-velocity profile, described above, 
that is consistent with the ranging solution. The bank angle command is 
generated from a vertical L/D command which is a function of a reference h/D, 
the difference between drag and drag reference, and the difference between alti- 
tude rate and altitude rate reference. The angle-of-attack profile is a func- 
tion of Earth-relative speed and consists of a series of linear and quadratic 
segments. The angle-of-attack profile is controlled through inputs in the mis- 
sion constants table. A complete derivation of all entry guidance equations can 
be found in reference 3* 




3.4.2 Entry Guidance Executive (EGEXbC) 

The EGEXEC function calls the other entry guidance functions in proper sequence 
and controls guidance phase transitions. In the normal mode EGFLG=0 or ^ 
guidance is divi ied into five phases: preentry, temperature control, equixibrium 

glide, constant and transition. Each phase, except preentry, computes a 

reference drag ac eration profile that is based upon the ranging requirement 

and the vehicle coi traints. 
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The computations performed during each guidance cycle vary with the guidance 
phase. The active guidance phase is defined by the integer variable ISLECT ac- 
cording to the following table. 

ISLECT Guidance Phase 

1 Preentry 

2 Temperature control 

3 Equilibrium glide 

4 Constant drag 

5 Transition 

The entry guidance phase transfer logic is summarized in table 3*4.2. 

The preentry phase commands a preselected bank angle in the EGPEP function. The 
preentry phase (ISELECT = 1) can be terminated by any of three conditions. 

a. Normal termination: The preentry phase is normally terminated at 5.66 

ft/sec2 by setting the ISLECT flag to 2 when the total acceleration exceeds 
the threshold value of ASTART ft/sec^. 

b. Alternate termination: If, at the threshold load factor level, the current 

constant drag level to reach the target is greater than the desired constant 
drag level (ALFM), ISLECT is set to 4. This could occur for an extremely 
short range case. 

c. Threshold load t actor level: If the current relative velocity is less than 

VTRAN, ISLECT is set to 5. 

The temperature control phase (ISLECT = 2) computes the required drag velocity 
profile during the high- temperature region of entry. The functions EGRP eind 
EGREF contain the detailed temperature control equations used by the entry guid- 
ance. During the temperature control phase, the drag-velocity reference 
trajectory is divided into two quadratic segments. Function EGREF determines 
which of these quadratics is to be used. Figure 3*4. 2-2 presents example qua- 
dratics for this phase. The quadratic for the higher velocity region is used 
when 

VE > VA1, and 

DRF = (DREF(2) - DREF(1)) (DHEF(2) - DHEF(1)) + (HDTRF(I) - HDTRF(2) GS1) > 0 

where DREF(l), HDTRF(I), DHEF(2), and HDTRF(2) are all computed in EGREF. If 
either test is failed, then the reference parameters RDTREF and DREFP computed 
for the lower velocity drag profile are used. The quadratic segment switching 
point is controlled by the mission load parameter GS1 . The quadratic switchirj 
logic, in the EGREF function, is similar to that used for guidance phase 
transitions . 
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The temperature control phase can be terminated in four ways: 

a. Normal termination: The temperature control phase normally transfers to the 

equilibrium glide phase, ISLECT = 3* Transfer is planned to occur when the 
drag reference profiles for the temperature control phase and the equilibri- 
um glide phase intersect, as illustrated in figure 3*4. 2-3. However, if the 
slopes of the two intersecting profiles are different, the bank-angle com- 
mand may Jump at the transfer point. Bank-smoothing logic provides a smooth 
roll-angle command history for minimum RCS fuel usage by transferring phases 
when the commanded bank angle is the same for both phases. This transition 
occurs before the intersection point of the drag profiles and is accom- 
plished by use of the following equation: 

If VE < VA, and 

DREFP < DREFP3 

then: 

ISLECT = 3 

where DREFP3 = DHEFP1 + GS2 (RDTREF - RDTRF1), and DREFP3 is confuted in the 
EGREF function. The value of GS2 determines the transfer point. 

b. Alternate termination: The guidance transfers to phase 3 when VE < VB1 if 

the bank- smoothing logic has not been satisfied by that time. Thus, if 

VE < VBl, then ISLECT = 3- 

c. Alternate termination for short ranges: For very short range targets, the 

desired constant drag level may be reached before equilibruim glide phase 
initiation. In this case, control is transferred directly from the tempera- 
ture control phase to the constant drag phase 

(ISLECT = l4) when 

VE < VCG + DELV, and 

DREFP > DREFPil 

where DREFP4 = T2 + GS3 (RDTREF + 2 HS T2/VE), and DREFP4 is computed in 
EGREF. Bank commano smoothing is provided through the constant GS3» 

d. Extremely short range termination: Transition to phase 4 occurs 

If T2 > ALFM 
then: 

ISLECT = 4 

wher^ T2 is computed in EGCOi'lN. 
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GS3. Transfer occurs when 
VE < VCG + DELV, and 
DREFP > DREFE'4 

where DREFP4 = T2 + GS3 (RDTREF + 2 HS T2/VE) 
and is again computed in EGREF. 

h Alterrate termination: For very long-range trajectories, the 

locitv at the intersection of the equilibrium glide and constant drag phases 
ifless than tift^ansition phase initiation velocity, VTRAN. When this 
occurs, the equilibrium glide phase transfers directly into the transition 
phase, ISLECT = 5. The transition to phase 5 occurs when 

VE < VCG + DELV, and 

VCG < VTRAN, and 

DREFP > DREFP5 

Where DKEFP5 = DF e (EFF - EEFh) Cl r 0S4 (BDIBEF - "OTBFI) la oo^^d in 
EOCOMN The variables EEF, C1, and RDTBFT are also computed In EGCOm. As 
lf“h;r phases. banlc-smooLlng logic is provided b, the msslon load eon- 

stant GS4. 

c. Alternate termination: For extremely short ’'f 6®®’ 

control phase, transition to the constant drag phase (ISLECT - 4) occurs ii 

T2 > ALFM 

where T2 is again computed in EGCOMN. 

The constant drag phase (ISLECT = 4) shapes the entry profile ® 

draff velocity profile to maximize the control system margins, 
coSalns the constant drag range prediction and controller 

stant drag phase terminates and transfers to the transition phase (^^LECT - 5) 
at a predetermined time before the transition energy level, ETRAN. Transfer 

occurs when 
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VE < VTRAN ♦ DELV, and 
DREFP > DREFP5 

where DREFP5 s DREFP5 in the equilibrium glide phase and the computed in EGCOMM. 

The transition phase is the final entry phase and is used to steer the Orbiter 
to the proper TAEM interface conditions. The transition phase and entry guid- 
ance are terminated at the start of the TAEM major mode. The entry- to-TAEM tran- 
sition logic is defined by 

If (ISLECT d 1 and VE < V_TAEM) then EG__END = 1 

To execute the entry guidance confutations properly, the following functions 
must be called in the sequence shown: 

a. The EGSCAlEHT function is called. 

b. On the first entry guidance pass (START 0), EGINIT is then executed. 

c. The EGCOMN function is called. 

d. The phase transition logic for the preentry phase (ISLECT = 1) and the first 
alternation termination test for the temperature control phase (if VE < VB1 , 
then ISLECT = 3) are executed within EGEXEC. The tests to transfer to phase 
4 if T2 > ALFM are also made at this time. 

The functions next called depend '>n the value of ISLECT: 

If ISLECr = 1, then perform EGPEP 

If ISLECT = 2 or 3, then perform EGRP, and then EGREF 
If ISLECT = 4, then perform EGREF4 
If ISLECT = 5, then perform EGTRAN 

After these functions have been executed, the remainder of the phase transition 
logic may be performed within EGEXEC at any time. The output commands are then 
computed by calling EGALPCMD ( angle-of-attack comnand) and the sequence 
EQGNSLCT, EGLODVCMD, and EGROLCMD (bank-angle command). If ISLECT s 1, 

EGGNSLCT and EGLODVCMD are bypassed; and only EGALPCMD and EGROLCMD are called. 
For all values of ISLECT, either the angle-of-attack or bank-angle command may 
be computed first. The entry guidance execution sequence is summarized in fig- 
ure 3. 4. 2-4. 

In the ETG "canned” mode EGFuG = 2, ISLECT is never allowed to be greater than 2. 
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3.4.3 Entry Guidance Scale Height (EGSCALEHT) 


The guidance function EGSCALEHT generates an altitude scale height (of 
atmospheric density), modeled on the 1962 standard atmosphere. This parameter is 
used in calculating the altitude rate reference term 


2D D cd 

RDTR.' = -HS ~ + -- 

VE D cd 


where D s drag acceleration 

D = time derivative of D 
cd = drag coefficient 
cd = time derivative of cd 

Empirical curve fits of the altitude scale height, HS, as a function of relative 
velocity, Ve, have been implemented into the entry guidance. 


3.4.4 Entry Guidance Initialization (EGINIT) 


The guidance function 
ance. In this routine 
time are computed. 


EGINIT serves as the initialization routine for entry guid' 
initial values are set and parameters calculated only one 


3 . 4.5 Entry Guidance Common Computation (EGCOW) 


The entry guidance contains several parameters used continuously throughout the 
guidance program. These are computed in BGCOMN and are such parameters as 
energy, EEF, the constant drag level to reach the target (T2), and the rate of 
change of T2, T2D0T. 


3 . 4.6 Entry Guidance Preentry Phase (EGPEP) 

In the Orbiter avionics system, the purpose of EGPEP is to generate a vertical 
L/D conmand (LODV) by means of the ILOAD parameter PREBNK. However, in 
for the MCC to simulate either an automatic or manual preentry phase, the LODV 
equation in the Orbiter avionics system should be replaced by the ba^-angle 
Input in hM 304 <j) by the flight controller. Also, the angle-of-attack com^nd 

issued in EGALPCMD should be overridden by MM 304a- The preentry phase, ISLECT 
= 1, is terminated by EGEXEC at a sensed total acceleration level equal to 
ASTART (currently 5.66 ft/sec^). 
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3 . 4.7 Entry Guidance Range Prediction (EGhP) 

THP FGRP function serves as the range predictor during the temperature control 
Z eoSinSriuHlSe phases. The range prediction is then used to determine 
the proper drag level during phases 2 and 3 to achieve the desired range at the 
entry-TAEM interface is only called when ISLECT is equal to 2 or 3 . 

In order to determine the proper drag-velocity profile for range control, a 
range prediction is made of the entire entry trajectory. This is accomplished 
b 7 using five drag-velocity segments; two during the temperature control phase, 
one during the equilibrium glide phase, one during the constant drag phase, and 
a constant range value during the transition phase, ^he temperature cont^ 
eauiiibnum glide, and constant drag phase range segments are computed in EGRP. 
The range value for transition, RPT, is computed in EGINIT and is constant in 
order tt provide a nominal transition range at the transition interface, VTRAN. 
Once the transition phase is entered, the transition range prediction is 
mouified in the EGTRAN function to meet the range requirements. 

The drag-velocity segments during the temperature control, equilibrium glide, 
and constant drag phases are anchored at specific velocity points as illustrated 

in figure 3 **<*T- 1 * 

The ranre for the temperature control phase is predicted along two quadratic 
drag-velocity segments anchored at VB1 and at VA1 . For the equilibrium g^-i^e 
L^nt, the ranic is predicted between VB1 and VCO. the oohputed interaection 
point between the equilibrium glide phase and constant drag phase, and for the 
constant drag segment, the range is predicted between VCG and VQ. In all cases, 
the entire drag-velocity profile is anchored at a drag level of 1^23 at a vel^- 
itv of VB Duriig the temperature control phase, VB is defined as VB1 ; and dur- 
ing the equilibrium glide phase, VB Is defined as the current relative velocity. 

There I’o re , 


if VE > VB1 


VB = VB1 


if VE < VB1 


VB = VE 


In the equilibrium glide phase, as the drag profile approaches the desired con- 
stant drag level, the locus of the equilibrium glide drag reference parameter 
may wander from a precise equilibrium glide profile shape trying to drive T2 to 
precisely ALFM, the desired constant drag level. In order to provide a more 
uniform equilibrium glide drag profile at the function point between the equilib- 
rium glide and the constant drag phases, the equilibrium glide reference 
is "frozen" when the rate of change of T2 is near zero. This is illustrated in 
figure 3 4.7-2. This is accomplished by freezing VB at the current value of VB 
when T 2 D 0 T becomes less than DT2MIN and when VE is less than VCG ♦ DELV. 


If Er.FLG is equal to 2, ETG canned mode, the ranging iteration of D23 is by- 
passed and Dc'3 is set equal to D23C in order to provide a "canned" drag veloc- 
ity profile for the ETG mode. 
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3.4.8 Rntry Guidance Reference Parameters (EGREP) 

ora.r to contnol tJ^OrtJLUr tTc^na 

parameters. The ° dependent part of the L/D reference parameter 

r j::‘te:;^*rn-'c:nirrs.r,:iuori!:. *uo. ph.». 1*1. fo.otio„ 1, 

only called when ISLBCT is equal to 2 or 3- 


Entry Guidance Constan t Drag Phase (EGRE^ 


3 . 4.9 — . 

The purpose of the depeLent^part of* the*L/D*referenoe for the 

«03trnrdr.IXo^: “Sirru^ticn ^ cau«. «1, -h» ISLECI 1. «.oal to M. 


Entiy Guidance Transition Ph ase (EGThAN) 


3.4.10 

The transition the^^constant^drag^phase^ transition phase 

reference level at the end of °°""„^”^,no«^Stes the correct drag-energy 

rofU^^^o r/r^r;nge ®-r 

r^?Sen^e*"?St; fincU^n is called only .*.en ISLECT is equal to 5. 


Entry Guidance Angle~Of ~Attack Function ^EGALPCMD). 


ontry .oloolt, regime finotton of relatlo. ..loo- 

of attack in each oegoeot is deflMd by a qmdratw lo^ 

tty Figure 3.4.11-1 shows a typical ^ 

iraies tSe rlexlblllt, avalUbl. In the profile sel«ltion. 

3. 4. ,2 Cn.r. Cutdanoe Cain Selec tion Funotlon (EOCIiaCT). 

THE eoONacr ration the gains 

C;i» if "braot^rdii: ^rit^tlon le.el and the dlffe«no. bet-een 
drag and drag reference. 
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3.4.13 Entry Guidance Lateral Logic and Vertical L/D Command Function 
(EGLODVCMD ) 

The purpoae of the EGLODVCMD function is to: 

Con 4 >ute the L/D reference parameter, ALDREP 

Compute a R feedback term to correct drag error biases caused by poor navi- 
gation (DLRDOT) 

Compute the vertical L/D connand from the controller equation: 

LODX = ALDREF + C16 (DRAG - DREFP) + C17 (RDTRF + DLRDOT - ROOT) 

Perform a velocity check to see if angle-of-attack modulation should begin 
in order to keep drag on the drag reference profile. 

Compute the baidc angle limit LMN and finally to compute the bank direction, 
RK2R0L. 



i 

! 

j 

\ 

i 

I 

! 


3.4.14 Entry Guidance Bank Command Function (EGROLCMD) 

The purpose of the EGROLCMD function is to generate a baink coomand for the auto- 
pilot and a bank reference parameter for display. The bank command is computed 
from the vertical L/D conmand parameter and the bank reference is computed from 
L/D reference. If the angle-of-attack modulation flag, ICT, is equal to one, a 
bank command bias is computed as a function of the ALPHA difference with respect 
to the nominal ALPHA schedule. 


3-4.15 Entry Guidance Data Flow Sunmary 

The data flow charts (figs. 3-4.15-1 and 3-4.15-2) present the data flow of all 
computea and stored parameters within the entry guidance. 


3.5 ENTRY AUTOPILOT FORMULATION 


3-5.1 Requirements Overview 

A simple 3-degrees-of-freedom autopilot is required to execute the bank and 
cingle-of-attack commands generated by the entry guidance. The autopilot require- 
ments are the same for both of the entry guidance modes. 

The autopilot, by means of a simple phase plane, generates the actual bank angle i i 

and angle of attack based on a maximum acceleration and rate limit in the bank ’ I 

and pitch channels, ihe sideslip angle ($) is always assumed to zero. i 

: ! 

The entry guidance generates a bank and angle-of-attack command necessary to con- ] 

trol the trajectory. The autopilot generates the Orbiter attitude response to J 


; } 
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these ""th^°a?Wudf reepinM S«t?irtS?cs!* tto'aSw^l” 

quency dynamics. Based on ito achieved during the next computer 

M wecirwlthln the “^uter cycle . ^ ”^ute 

ih:TrafeTtc “'dySSr. .™i fcr the n«tt lnte«r.tich step dur- 

ing entry. 


3 . 5.2 Autopilot Executive (DAP3P1 


The autopilot ex^uUve rcuUce PAP^ Tp l^er^torPHSPUh 

r SlunuSTcnrriranrUrerru^rconthcl. t^apectlvely. 


DAP 3 D also computes the “utlter attit^e with respect to '|ja JS^eiera- 

r'o?ihr^^ortoT„irt?iin^r^toaS^^^^^ 

^■d'tobtoTreTreSrtfthe totor^l^parameter definitions. Appendix B pres- 
:"?s toi fc;^la?ton ftow SLrts, and Appendix D presents the IBM structured 

flow charts for the autopilot. 


3.5*3 Autopilot Phase Plane (PHSPLNj^ 


The PHSPLN routine is a simplified phase Plane ^^ai i® us^ to 

r ±'bbm"c: r^ansr aXudrd— rorrPAr^ 

to be selected by the user. 


3.5.4 Autopilot Data Flow Summary 

The data flc charts (figs. 3.5.»-' snd 3 - 5 . 4 - 2 ) present the daU flow of all 
computed and stored parameters with the autopilot. 


3.6 TARGETING ROUTINE (EGKT) 

SrXt?SSSS3S=lS^^ 
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point is converted into an Earth-fixed position vector, and the great circle 
range to target is computed between the vehicle and this tangent point. The arc 
length is then computed from the tangent point around the alinement cone to WP1, 
the straight in approach point. The range to target is then computed as the sum 
of the great circle range to the tangent point, the arc length around the 
alinement cone, and the distance between WP1 and the runway threshold point. The 
azimuth error is computed as the dilferenoe between the vehicle Earth-relative 
azimuth and the heading to the heading alinement tangent point. Change 8 incor- 
porates the required modifications to the targeting logic for optional TAEM 
targeting (OTT) capability. This capability allows selection of either an 
overhead (far HAC initially) or straight-in (nearest HAC) approach for both pri- 
mary and secondary runways independently by MED input. The desired mode is 
obtained by setting the input flag OVHD to one or zero to to obtain either an 
overhead or straight-in approach, respectively. This flag is a two-element 
array wherv the first element denotes primary runway and second element denotes 
secondary runway (secondary is required only for midtargeting and redesignati n 
options). The selected runway designation flag (RWID) is required as input f r 
proper initialization at time of redesignation. RWID = 1 if the primary runway 
is selected and RWID = 2 if a redesignation to the secondary runway has 
occurred. The velocity to switch HAC targeting overhead status is nominally 
zero but may be input via MED to simulate a manual HAC toggle during the entry 
trajectory. The midpoint targeting is needed to minimize the delta range anc az- 
imuth error for a redesignation due to low altitude winds. Midpoint targeting 
is accomplished by targeting to both primary and secondary runways and averaging 
the range and azimuth error values and is terminated when 

a. The crew selects the desired runway 

b. The relative velocity satisfies a preset velocity limit and the primary 
runway is selected by default 

The flow charts for the targeting routine are found in appendix C and are 
subdivided into the following functions: 

f 

. 3.6.1 EGRT-E.XEC, Targeting Executive Logic 

This routine performs the necessary initialization of targeting parameters and 
calls the EGRT sequence controller. If midpoint targeting is active, the EGRT 
sequence is recycled for the secondary runway, and the resulting ranges and azi- 
muth errors are averaged. The EGRT sequence controller executes the following 
sub functions for the specified input argument list. 

3.6.2 EGRT-dHACRC, Center of Heading Alinement Cone - Runway Coordinates 

This routine computfS the position of tfie center of the heading alinement cone 
in runway coor'd inates . 
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3 . 6.3 EGBT-CHACRC. Center of Heading Alineae nt Cone in Earth-Fixed Coordinates 

This function transfoms the center of the heading alinement cone vector to the 
Earth-fixed coordinate system. 


3.6.M EGRT-BV. Bearing of the Vehicle 

This function computes the bearing of the vehicle based on the current vehicle 
Earth-fixed position vector. 

3.6.5 EGRT-BVCHAC. Bearing to Center of th e Alinement Cone 

This function computes the bearing from the vehicle to the center of the heading 
alinement cone. 

3.6.6 EGRT-COSTHETA. Great Circle Arc 

This function computes the great circle arc between the vehicle and the center 
of the heading alinement cone. 


3.6.7 EGRT-DWP 1 . Distance to WP 1 

This function computes the range to the tangency point on the heading alinement 
cone . 


3.6.8 EGRT-DVNEP. Range- to-Threshold Point 


This function computes the heading 
alinement cone, the distance around 
the runway threshold point. 


to the tangency point on the heading 
the alinement cone, and the final range to 


3.6.9 EGRT-DELAZ. Azimuth Error 

This function computes the azimuth error between the vehicle heading and the 
heading to the tangency point of the alinement cone. 
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TABLE 3.3-1.- ENTRY GUIDANCE INPUT DATA 
(a) Input parameters 


Symbol 

Description 

Source 

Unit 

ALPHA 

Angle of attack 

Autopilot 

deg 

DELAZ 

Current azimuth error 

EGRT 

rad 

DRAG 

Current drag acceleration 

State vector 

ft/sec2 

EGFLG 

Guidance mode flag 

Entry processor 
executive 

n.d. 

HLS 

Altitude above runway 

State vector 

ft 

LOD 

Current lift/drag ratio 

Aerodynamics 

n.d. 

RDOT 

Current oblate Earth altitude rate 

State vector 

ft/sec 

ROLL 

Current bank angle 

Autopilot 

rat 

START 

Initialization flag 

Entry processor 
executive 

n.d. 

TRANGE 

Current range to runway 

EGRT 

n. ml. 

VE 

Current relative velocity 

State vector 

ft/ sec 

VI 

Current inertial velocity 

State vector 

ft/sec 

XLFAC 

Current load factor 

Aerodynamics 

ft/sec^ 

MM304(|) 

Preentry bank angle 

Med input 

deg 

MM304a 

Preentry angle of attack 

Med input 

deg 

MEP 

Left-hand HAC select flag 

n.d. 
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TABLE 3.3-1.- ENTRY GUIDANCE INPUT D.\TA - Continueo 
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TABLE 3.3-1.- ENTRY GUIDANCE INPUT DATA - Continued 
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TABLE 3.3-1.- ENTRY GUIDANCE INPUT DATA - Continued 



CT17M2 Value of CT17MN when ICT = 1 0.00133 sec/ft 


TABLE 3.3-1.- ENTRY GUIDANCE INPUT DATA - Continued 




TABLE 3.3-1.- ENTRY GUIDANCE INPUT DATA - Continued 
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TABLE 3-3-2.- ENTRY GUIDANCE OUTPUTS 


Symbol 

Description 

Unit 

Destination 

ALPCHD 

Angle-of-attack oonmand 

rad 

Autopilot 

ROLLC 

Bank conmand 

rad 

Autopilot 

drefp 

Drag reference 

ft/sec2 

Display 

DRAG 

Actual drag 

ft/sec^ 

Display 

ROLREF 

Bank reference 

deg 

Display 

ISLBCT 

Guidance phase Indicator 
relative 

n.d. 

Display 

VCG 

Velocity at start of constant 
drag phase 

ft/ sec 

Display 

VRR 

Velocity at first bank 
reversal 

ft/ sec 

Display 

BOWD 

Energy overweight 

ft 

Display 

BBI 

Entry evaluation indicator 

n.d. 

Display 

RC176G 

First roll command after 0.176g 

deg 

Display 
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TABLE 3.3-3.- ENTRY GUIDANCE INTERNAL PARAMETER DEFINITIONS 




Symbol 

Description 



ACLAM 

Maximum allowable alpha 



ACLIM 

Minimum allowable alpha 



ACMD1 

Scheduled angle of attack 



ALDCO 

Temporary variable In phase 3 reference parameters 

ALDREF 

Vertical L/D reference 



ALPCMD 

Angle-of-attack command 



ALPDOT 

Rate of change of ALPCMD 



ARG(I) 

Cosine of commanded bank angle 



ARG(2) 

Cosine of unlimited band command 



ARG(3) 

Cosine of bank reference angle 



A2 

Temporary variable used in computing 

range 

and updating D23 

CAG 

Pseudoenergy/mass used in transition 

(L/D) 

reference 

CQUl ,2) 

Constants in Ith temperature control D-V quadratic 

CQ2(1,2) 

VE coefficients in Ith temperature control 

D-V quadratic 

CQ3(1 ,2) 

VE^’ coefficients in Ith temperature 

control D-V quadratic 

Cl 

dD/dE in transition 



C16 

d(L/D)/dD 



C17 

d(L/D)/dH 



C2 

Component of L/D reference 



C4 

Reference altitude rate term 



C20 

d ALPHA/dCd gain 



DD 

Drag - DREFP 



DDS 

Limited value of l)D 
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Symbol 

IO)P 

PELAM*' 

DEL ALP 

DLRD07 

DLIM 

DLZRL 

D»DD 

DREP'(1,2) 

DRF5P 

DREFPT 

DREFP1 

DREFP3 

DREFPii 

DREP'?5 

DRF 

DX(1,2) 

DZOLD 

DZSGN 

D231 

KEF 

!-7)7RF ( 1 ) 

I ALP 


TABLE 3*3**3«~ Continued 


Description 


Past value of OD 
Delta ALPHA from scheclule 
Command ALPHA increment 
R feedba(dc term 

Maximum value of DREFP in transiticm 

Test variable in bank angle computation 

Derivative of range wrt drag 

DREFP for Ith ten^ierature control D-V quadratic 

Drag reference used in controller 

DREFP-DF in transition phase 

DREFP in eouilibrium glide 

DREFP test value for transition to phase 3 

DREFP test value for transition to phase 4 

DREFP test value for transition to phase 5 

Test value for transition to D23-VB1 quadratic reference 
parameters 

Normalized values of DREFP 
Previous value of DELAZ 
Change in DELAZ 
First updated value of D23 
Energy /mass 

Intermediate calculation of temperature control R ref 
ALPCMD segment counter 

I 
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TABLE 3 . 3 - 3 .- Continued 


Symbol 

Description 

ICT 

Alpha modulation flag 

ISLECP 

Past value of ISLECT 

ISLECT 

Entry guidance subphase counter 

ITRAN 

Transition initialization flag 

LMFLG 

Saturated roll command flag 

LMN 

Maximum value of LODV 

LODV 

Vertical L/D command 

LODX 

Unlimited vertical L/D command 

Q(i,2,3) 

DREFP/VE in temp control phase 

RCG 

Constant drag phase range 

RCG1 

Constant component of RCG 

RDEALF 

Roll bias for alpha modulation 

RDTREF 

Altitude rate reference 

RDTRF1 

RDTREF in Phase 3 

REQ1 

Equilibrium glide range x D23 

RER1 

Transition phase range 

RF(1,2) 

Ith range segment in temp control phase x D23 

RFF1 

Temperature control range x D23 

RDTRF 

Altitude rate reference corrected for Cd 

ROLLC(I) 

Roll angle command about body axis 

R0LLC(2) 

Unlimited roll command 

R0LLC(3) 

Rolref 

RPT 

Desired range at VQ 
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TABLE 3>3“3-- Concluded 


Symbol 


R231 

START 

T1 

T2 

T2DOT 

T20LD 

V(1,2,3) 

VB2 

VCG 

VE2 

VF(1,2) 

VSAT2 

VTRB 

VX(1 ,2) 

XLOD 

YL 

ZK 


Description 

Phase 2 and 3 range x D23 

Entry guidance first pass flag 

Equilibrium glide vertical lift acceleration 

Constant drag level to reach target 

Rate of change of T2 

Old T2 value 

Velocity sampling points for temp control numerical range 
prediction 

Vb2 

Phase 3-^ boundary velocity 
VE2 

Upper velocity bounds for Ith temperature control range 
segment 

VSAT^ 

R feedback velocity lockout 

Velocities where dD/dV=0 in Ith temp control D-V quadratic 
Limited left/drag ratio 
Maximum heading error absolute value 
R feedback gain 


J 
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(a) Input parameters 
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m 
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Symbol 

Description 

Source 

Unit 

ALPWWD 

Angle of attack with winds 

Integrator 

rad 

ALPHAP 

Pr6vious pass angls of attack 

Entry integrator 

rad 

ALPCMD 

Angle-of-attack command 

Entry guidance 

rad 

PQR(IVI) 

Initial bank rate 

Entry integrator 

rad/ sec 

PQR(IV+2) 

Initial alpha rate 

Entry integrator 

rad/ sec 

ROLLC 

Bank command 

Entry guidance 

rad 

ROLLP 

Previous pass bank angle 

Entry integrator 

rad 

RRPAST 

Previous pass bank rate 

Entry integrator 

rad /sec 

PRPAST 

Previous pass angle-of-attack rate 

Entry integrator 

rad /sec 

EGFLG 

Guidance mode flag 

Executive 

n.d. 

F5M304(j, 

Preentry bank command 

Med input 

rad 

MM304a 

Preentry angle-of-attack command 

Med input 

rad 

ROLL 

Current ba.^k angle 

Integrator 

rad 

DELTT 

Integration time step 

Executive 

sec 
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TABLE 3.3-4.- AUTOPILOT INPUT DATA 


(b) Input constants 


Symbol 

Description 

Value 

Unit 

Class 

DTK 

Degrees to radians 

0.017453292 

rad 

f 

KONTROL 

Initialization flag 

15 

n.d. 


PA 

Maximum angle-of-attack 
acceleration 

0.0174533 

rad/sec^ 

f 

PADB 

Angle-of-attack attitude 
deadband 

0.00174533 

rad 

f 

PRM 

Maximum angle-of-attack 
rate 

0.0872665 

rad/ sec 

f 

RA 

Maximum bank acceleration 

0.029671 

rad/sec^ 

f 

RADB 

Bank attitude deadband 

0.043633 

rad 

f 

RADB2 

Bank attitude deadband for 
EGFLG=2 

0.00174533 

rad 

f 

RAD 160 

160 degrees converted to 
radians 

2.7925268 

rad 

f 

RAD180 

180 degrees converted to 
radians 

3.1415927 

rad 

f 

RAD 360 

360 degrees converted to 
radians 

6.2831853 

rad 

f 

RRM 

Maximum bank rate 

0.0872665 

rad/ sec 

f 


fl 
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TABLE 3.3-5.- AUTOPILOT OUTPUTS 


Symbol 

Description 

Unit 

Destination 

ALPWND 

Updated angle of attack with 

rad 

Integrator 


winds 



BANK 

Updated bank angle 

rad 

Integrator 

PR 

Updated angle-of-attack rate 

rad /sec 

Integrator 

RR 

Updated bank rate , 

rad/ sec 

Integrator 
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TABLE 3.3-6.- AUTOPILOT INTERNAL PARAMETER DEFINITIONS 


Description 

Bank about the velocity vector 
Difference between ROLLC and ROLL 
Direction of acceleration 
DAP 3D time step 

Average value of EFRATl and EFRAT2 
DIFF rate 

Past value of EFRATl 
Cycle frequency of DAP3D 
First and third pass flag 
Flag to establish constant set 
PHSPLN initialization flag 
Current bank rate 
Current bank angle 


\ 

i 

1 

j 

i 

i 

i 

1 


j ' 
i ' 



Symbol 

BANK 

DIFF 

DIR 

DTIM 

EFRATE 

EFRATl 

EFRAT2 

ICPPLN 

ICTTRN 

INTRY 

KONTRL 

RR1 

ROLL 
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TABLE 3.3-7.- TARGETING ROUTINE INPUT DATA 


(a) Input parameters 


Symbol 

Description 

Source 

Unit 

AZRH 

Bearing from true north of runway 
* X axis 

Landing site table 

rad 

XYZE 

Vehicle position vector 

Integrator 

ft 

XYZED 

Vehicle velocity vector 

Integrator 

ft/sec 

SRAZ 

Secondary runway azimuth 

Landing site table 

rad 

M 

Vehicle mass 

AERO 

slugs 

(oVHD)a 

Runway approach mode flags 

MED 

n.d. 

RWID 

Selected runway ID flag 
(pri = 1, sec = 2) 

site selection 

n.d. 

(rec)*^ 

Greenwich to runway transformation 
matrices 

site selection 

n.d. 

(rLS)c 

Runway position vectors 
(Greenwich coord.) 

site selection 

ft 

IFP 

First pass flag (= 0 initially) 

EXEC 

— 

VTOGL'* 

Velocity to toggle OVHD/STRT HAC 
status 

MED 

ft/ sec 


*MED input to select approach mode (overhead = 1, straight in = 0) (dimensioned 
(2) >*iere subscript is runway ID flag). Initialized to 1 in mission constants 
table . 

^Dimensioned (3x3x2) where last index is runway ID. 
dimensioned (3 x 2) where last index is runway ID. 

MED input to simulate manual HAC toggle. Initialized to zero in mission con- 
stants table. 
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TABLE 3 •3-7.- Concluded 
(b) Input constants 


Symbol 

Description 

Value 

Unit 

Class 

RX22 

Polar radius^/equatorial 

radius^ 

0.9933065782 

n.d. 

f 

VMIDPT 

Velocity limit for midpoint 
targeting 

TBD 

fps 

m 

XHACL 

Low mass RW to HAC distance 

-352A5. 

ft 

m 

XHACH 

High mass RW to HAC distance 

-35245. 

ft 

m 

WTGS1 

Mass threshold value 

8000. 

slugs 

m 

HI 

Initial HAC radius 

20000. 

ft 

f 

PSHARS 

Initial HAC turn angle 

270. 

deg 

f 

RFO 

HAC radius on final 

14000. 

ft 

f 

R1 

Linear term in spiral 
radius eqn. 

0. 

ft/deg 

f 

R2 

Quadratic terra in spiral 
radius eqn. 

.093 

ft/ deg^ 

f 

A3T0L 

HAC turn angle tolerance 
value 

-.003 

deg 

f 

DTR 

Deg to rad conversion const. 

.01745329 

rad /deg 

f 

RTD 

Rad to deg conversion const . 

57.29578 

deg/ rad 

f 
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Symbol 

TRANCE 

DEUZ 

RCHMAG 

PSHAT 

RTURNT 

YSGNT 


TABLE 3.3-8.* TARGETING ROUTINE OUTPUTS 


Description 

Unit 

Destination 

Range to target 

n. ml. 

Entry guidance 

Azimuth error 

rad 

Entry guidance 

Radius of landing site 

ft 

Deorblt processor 

HAC turn angle 

deg 

TAEM guidance 

HAC radius 

ft 

TAEM guidance 

R/L HAC indicator 

n.d. 

TAEM guidance 
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TABLE 3 . 3 - 9 .- TARGETING ROUTINE INTERNAL PARAMETER DEFINITIONS 


Symbol 


Description 


RLS 

RC 

HACEF 

BARCC 

CTHVC 

BAKWPI 

DARC 

PSI 

REC 


Landing site in Earth-fixed coordinates 

Center of heading alinement circle in runway coordinates 

Center of heading alinement circle in Earth-fixed coordinates 

Heading to center of alinement circle 

Cosine (.angle between vehicle and HAC) 

Heading to tangent point on HAC 
Distance around heading alinement cone 
Heading of vehicle 
Runway to Earth-fixed matrix 
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TABLE 3.W.- GUIBAJjCE PHASE SELtCIIOh LOGIC 



'1 




Name: Greenwich true of date (geographic). 

Origin: The center of the Earth. 

Orientation: The X^-Yg. plane is the Earth's true*of-date equator. 

The Zg-axia is directed along the Earth's true-of-date 

rotational axis and is positive north. 

The '•‘X -axis is directed toward the prime meridian. 

E 

The Y -exie complete* * right-handed *y*Cem. 

E 

Charecteriatics: Rotating, right-handed, Carteaian. Velocity vector* 

expreaaed in thia ayatem are relative to a rotating reference 
frame fixed to the Earth, whose rotation rate* are expressed 
relative to the Aries-mean-of-l930 system. 


Figure 3.2-1.- Greenwich true of date (geographic). 
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NORTH 



Name: Runway coordinate system. 

Origin: Runway center at approach threshold. 

Characteristics: Rotating, Earth-referenced. 

Description: Zj^-axis is normal to the ellipsoid model through the runway 

centerline at the approach threshold and positive toward the 
center of the Earth. X^-axis is perpendicular to the Z,^-axis 
and lies in a plane containing the Zj^-axis and the runway 
centerline (positive in the direction of landing). 

Y„-axis completes the right-handed system. 

K 

is the runway azimuth, measured in the plane from 

true north to the +Xp-axis (positive clockwise). 


Figure 3.2-2.- Runway coordinates. 
40 
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Figure 3. 4. 2-3.- Bank-angle smoothing logic between guidance phases 
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Figure 3. 4. 7-1.- Drag-velocity segments used In range predictions 



VE 


Figure 3. 4. 7-2.- Freezing of equilibrium glide profile. 
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Figure 3.4.15-1.- Entry guidance external data flow suitmary. 


47 






r 




48 


.44 






STORAGE 


EGPEP 



Figure 3.4.15-2.- Continued. 


49 









STORAGE 


DREFP, RDTREF. C2 
ITRAN 


ALPCMD, ALPDOT <4 


C16, C17< 


OLROOT, ALPHA, 
ALDREF. LOOX, EE I 
LODV, RK2ROL, 
DELALP, ICT -4 


DRAG, XLFAC 
TRANGE, VE 


JTRAN, DREFP, |EEF,XLOD 


HS. VE2, CAG 


VE, MM304°= 


DRAG, ICT 


DELAZ, VE, 
drag, RDOT 


ROLL 


RDTREF, ISLECT, VE2 

VSAT2, C2, HS, C16 
Cl 7. DREFP, XLOD 
RK2R0L, VSAT2 
ALPCMD, ALPDOT, T1 
ICT 


EGTRAN 


IF ISLECT - 5: - 
COMPUTE REFERENCE 
PARAMETERS 


EGALPCMD 


COMPUTE AN6LE-0F- 
ATTACK COMMAND 


EGGNSLCT 


IF ISLECT. > 1: 
COMPUTE L/D 
COMMAND GAINS 


EGLODVCMD 


IF ISLECT>1: 

COMPUTE L/D, 

ROLL DIRECTION. AND 
ALTITUDE RATE COR- 
RECTION FOR Cq 


Figure 3.4.15-2.- Continued. 
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Figure 3.4.15-2.- Concluded. 
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Figure 3. 5. 4-1.- Autopilot external dau flow. 
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Figure 3. 5. 4-2.- Autopilot internal data flow. 
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APPENDIX A: 


ENTRY GUIDANCE FLOW CHARTS 


The following flow charts define the entry guidance formulations. 
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Figure A-1.- EGEXEC, entry guidance executive. Page ' ^ 
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Figure A-1.- Continued. 
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Figure A-1.- Continued. ^*9® 3 of 4 
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Figure A-1,- Concluded. 
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Figure A-3.- E6INIT, initialization. 
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Figure A-3.> Concluded. 
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Figure A-6.- EGRP, range prediction. Page 1 of 4 
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Figure A-6.- Continued. P*96 2 of 4 
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Figure A-6.- Continued. 
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Figure A-9.- EGTRAN, transition phase. 


Then »*♦♦♦*♦♦•♦*♦♦♦**• 

« OREFP = OLIM; ♦ 





Figure A-9.- Concluded. 
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Figure A-10.- E6ALPCMD, angle-of-attack command. 









Figure A-11.- EGGNSLCT, gain select. 
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I' COMPUTE VEATICAl L/D COMMAND (LODV) 

t 

f 

■ A4i4 = EXP(.(VE - CDDOTl) i CDDOT?); • 

. • • 

^ • COCAL = CDDOT4 ♦ ALPCMD (CDDOT5 ♦ CDDOTb • 

^ • • 

• ALPCMD) ♦ CDDOT i A 44; * 

• • 

• CDDOTC = CDDOT7 (DRAG ♦ GS ROOT / VE) • 

• • 

• A44 ♦ ALPDOT (CDOOTb ALPCMD ♦ CDDOH); • 

• • 

• C4 s HS CDDOTC / CDCAL; • 


i 

I 

I 

.1 


• 1CT s 1; • 

t t 

• C^O a MIDVAL(C21, C22 ♦ C23 VE» C24); • 


then ••*•*•*••§••••§•§*■•••••§••••••••••••§••• 

• IF VE < VC20 •— • C20 a MAX(SCALAB(C25 ♦ C2 d VE , C27)); • 


1 * DELALP - MIOVAL(CDCAL ((OREFP / DRAG) - • 

I * * 

I • 1.) / C20. OLAPLM, - DLAPLM); * 

then 

• IF ISLECT = 5 • • T1 a GS (VE2 / VSAT2 - 1.); • I 


I IF ABS (DRAG - DREFP) < DDMIN \ IDELALP * 0.0 I 

1 t \ I 




• IF VE < VNOALP •- 


-• IF DRAG > s DREFP OR VE < VALMOD OR ICT a 

• 1 


....... 


• AI.DREF a TW DREFP ♦ (2. RDTREF ♦ C2 • 1 

t • 1 

• HS) / VE; • I 

• f I 

• RDTRF a RDTREF ♦ C4; • I 

• * I 

• DD a drag - DREFP; • | 


I 

I 

i 

then 

• IF VE < VRDT • • DDS » MIDVALU'D, -DDLIM, DDLIM); 

• ZK a 2K1; 


Figure A- 12.- EGLODVCMD, lateral Icglc and verticil L/D coaunand. 
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Figure A-12.- Continued. 



THEN 






Figure A-12.- Concluded. 







COMPuffc KJtc CUNNANOS AND REfEkcNCL ANGLES 


IF(RK2R0L RK2RLP < 0 AND IVRR - 0^^ 


THEN 


VRR « VE 
IVRR • I 


0 AHG « LUUV / <LiiC; « 
* ^ « 

• ARG • LOOX / XLUU; * 

* ^ * 

♦ • 
« AKG « ALUREF / XLJOi • 
♦ J • 


• DO FOK I • i ru 3; •- 


IF AbS(A><u ) >* 1. 
I 


THEN <•*••••*•••••••••••••♦••• 

« « SIGNiARG XL2UI; * 

• I I * 


• RULLC « RKZROL ARCCOS(&HG t. ' t 

• I I ' ♦ 



IF AdSIROLLC ) > RLk AND VE < VEROtC • 

I • 


Then 



-• ROUC ■ «LN SICNIROLLC M J 


♦RGLREF FUR DISPLAY IS ROLLC(3) 


Ffqure A- 13.- E6R0LCM0, roll coimand. 
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appendix B: ENTRY AUTOPILOT FLOW CHARTS 

T,» fouowins n», Chart, Cflne t« cntr, autopilot fo»ulationa, 


Function 

DAP3D 

PHSP^JJ 


Figure 


Miimher of flow charts 
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Figure B-1.- UARJU. 
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Figure B-2.- Continued. 
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T 


E6RT-EXEC 


IF IFP . 0 


RWIDO • RWIO 


IF RMID t RHID^ 



IF V£ > VmOPT 


IF RWIO • I 


VE • |XYZED| 



IMID ■ 1 



THEN 

— 

RAZ . AZRH 

1 1 


Ll^Lj 

RAZ . SRAZ 


CALL EGRT-SEQ 
OUTPUT LIST- 

rang. DELAZ. YSGN(RWIO). RTURN(RWIO). PSHA(RWIO) 


IF IHID • I 


I 'ISONT • TSGN(RWIO) 

PSHAT . psha|r«i5{ 

RTURNT • RTURN(RHIO) 


CALL EGRT-SEQ * 

INPUT LIST; 0VH0{2). RLSfl ?) Bren 1 

RTURN(2), PSHA(2), IMIO ‘ . 

STRANG. SDELAZ. YSGN( 2 ). RTURN( 2 ). PSHA( 2 ) 


DELA 2 • (DELA 2 ♦ S 0 ELA 2 )/ 2 . 
IRANC • (TRANG ♦ STRANG),' 2 . 


fljur. C.I.. EGRT-EXEC, t.rj.tlnj 


•»»cutlv« lojtc. 











EGRT-SEQ 


INPUT LIST: OVHO, RLS, REC. RA2, YSGN. RTURN, PSHA, IMIO 
OUTPUT LIST: TRANG. OELAZ» YSGN. RTURN, PSHA 













Figure C-3.- E6RT-CHACRC, center heeding allnement 
cone - runway coordinate system. 
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E6RT-CHACEFC 




[HACEF] = 

[REC]’^[RC] 




Fi9ure C-4.- EGKT-CHACtFC , center he.idiny alinement circle - 

Earth- fixed coordinates. 
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E6RT-6V 



RCCEF(I) - RLS(I) + HACEF(I) 


vSl(2) * xmlalRSpIn ■ ^^^^p)RCCEF(2) 

v*«^(3l . xyzl{?)gg|^g{ : 

«VEHM(i = /xVZE(1)2 + XYZE(2)2 + X»2E(3)2 
T3 = VN0Hi'1(3)RVEHMG 

“ ^«F(l)2tRCCEF(2)2*RccEF(3)2 

SINB = RTURN/RCHf-lAG 


Figure C-5,- EGRT-BV, bearing of vehicle. 
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Figure C-6.- EGKT-BVCHAC, bearing to center of alinement circle. 
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C-7 i 





Figure C-7.- EoRT-CUSTHETA, great circle arc. 
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EGRT-DWPl 


STHVC =7mCRT01 - CTHVC^ 

CTVWPl » CTHVC + MCRT02 CTHVC SINB 

SBARCR * SINB/STHVC 

TEMP = CTVWPl SBARCR 

TEMP = AMIN1(1.,AMAX1(TEMP,-1.)) 

A2 » ACOS(TEMP) 

TEMP = CTVWPl 

TEMP = AMIN1(1..AMAX1(TEMP.-1.)) 
DVEWPl = ACOS(TEMP)RCHMAG 


Figure C-8«- EGKT-UWPl , distance to WPl. 
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C-9 




EGRT-DVNEP 




IF PSHA > PSHARS + 1 . 
OR 

A3 < A3T0L 
OR 

YSGN f SIGNY 


A3 = A3 + 2 PI 



Figure C-9.- EGRT-DVNEP, range- to- threshold point. 
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Figure C-10.- EGRT-DELAZ, azimuth error. 
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APPENDIX D: IBM AUTOPILOT FLOW CHARTS 

The flow charts presented in this appendix represents IBM's response to the re- 
quirements set forth in this document. The autopilot model for the MCC will be 
programed by IBM from the flow charts contained in this appendix. 





Figure D-1.- Entry autopilotv 







Figure 0-1.- Continued. 







Figure D-1.- Continued. 





Continued 
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Figure D-).. Continued. 






DIFF Is Shorter Than 



Continued. 

















Figure D-1.- Conclude^l. 



